Long-Term Transit Timing Monitoring and Refined Light Curve Parameters of 

HAT-P-13b 

Benjamin J. Fulton^, Avi Shporer^'^, Joshua N. Winn^, Matthew J. Holman'^, Andras Pal^'^, J. 

Zachary Gazak^ 



ABSTRACT 



We present 10 new transit hght curves of the transiting hot Jupiter HAT-P-13b, 
obtained during two observational seasons by three different telescopes. When combined 
with 12 previously published light curves, we have a sample consisting of 22 transit 
light curves, spanning 1,041 days across four observational seasons. We use t his sample 



to ex amine the recently observed large-amplitude transit timing variations (jPal et al. 



20111), and give refined system parameters. We find that the transit times are consistent 



with a linear ephemeris, with the exception of a single transit time, from UT 2009 
Nov 5, for which the measured mid transit time significantly deviates from our linear 
ephemeris. The nature of this deviation is not clear, and the rest of the data do not 
show any significant transit timing variation. 



Introduction 



HAT-P-13 (jBakos et al.l l2009l ) is among the brightest stars (^=10. 6 mag) hosting a multi- 
planet system containing a transit ing planet, HAT-P-13b (Mp^ = 0.85 Mj, Rp}j = 1.3 Rj, 



Bakos et al. 



200S 



Winn et al.ll20ld ). A transit depth of ~1% and a short orbital period of ~2.92 



days allow for many transits that can be observed by small ground-based telescop es, makiiig this 



system a good subject of observationa l studies related to planetary systems (e.g., iMardling 
Winn et aLlboiol : IPavne fc Fo^boilh . 
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HAT-P-13b was discovered as a transiting planet by iBakos et al.l (1200911 . who also identified 



a second planet in the system, HAT-P-13c (Mp^c sin v = 14.3 Mj^ 
in an eccentric orbit (e ~ 0.7), with a period of about 1.2 years 



Winn et al. 



Winn et al 



2010), moving 
(|20ld l gathered 



additional radial velocity (RV) measurements and identified a third low mass companion in the 
system, possibly a thir d planet, whose per iod is currently unknown but expected to be a few years 
or longer. In addition, I Winn et al.l (|20ld ) i dentified that HAT -P-13b orbit is likely to be aligned 
with t he host star's equator. As indicated by I Winn et al.l (120101 ) . based on the analysis of lMardling 
(j2010l ). a spin-orbit alignment of HAT-P-13b suggests a small mutual orbital inclination of planets 
b and c, suggesting in turn the possibility that planet c is also transiting. So far no transits of 
HAT-P-13c have been detec ted, although soin e attempts were made to look for it during the 2010 
predicted conjunction time (jSzabo et al.ll201Cl ). 



In a recent paper iPal et al.l (120111 ) analyzed transit timing of HAT-P-13b from four observa- 
tional seasons, 2007/2008 (hereafter Season 1), 2008/2009 (hereafter Season 2), 2009/2010 (here- 
after Season 3) and 2010/2011 (hereafter Season 4), and identified a deviation of the transit times 
from the predicted times during the last season, of about 0.015 day. In principle, this kind of long 
term transit timing variation (TTV) could be due to the presence of another planet in the system, 
in a large eccentric orbit, as described for example bv lAgol et al.l (l2005l. their Section 4 ). The 
known orbit of HAT-P-13c does not match the TTV patter n iden tified bv lPal et al.l (j201ll ). but it 
could be due to the third planet suggested by I Winn et al.l (|2010l ) , or a further companion in the 
system. 

We have set out here to study the suggested TTV signal in more detail. We present 10 new 
HAT-P-13b transit light curves, 5 from each of th e Seas ons 3 and 4, and combine them with the 
12 light curves which were available to iPal et al.l (j201ll ). Therefore, our analysis is based on a 
total of 22 transit light curves, either partial or complete, from four consecutive observational 
seasons, including a single light curve form Season 1 and 6-8 light curves per season for Seasons 
2-4. Overall, our data span 1,041 days. Obtaining the new data and photometric processing is 
described in Section 2, and in Section 3 we present our transit light curve analysis. We discuss our 
results in Section 4, and give a short summary in Section 5. 



2. Observations 

Our 10 new transit light curves of HAT-P-13b from Seasons 3 and 4 were obtained at three 
observatories. A brief description of the three telescopes and instruments used is given in the 
following paragraphs. 

Fred Lawrence Whipple Observatory (FLWO): FLWO is located on Mount Hopkins, near 
Amado, AZ. We used the KeplerCam 4kx4k Fairchild CCD486 mounted on the FLWO 1.2m 
telescope. KeplerCam has a pixel scale of 0'.'62 pixel"^ (2x2 binning), and a 23fl x 23fl field of 
view (FOV). All four new FLWO light curves, two from each of Seasons 3 and 4, were obtained in 



- 3 - 



the SDSS-z' filter. These observations were conducted with the telescope nominally in focus, but 
the optical characteristics of the telescope create a relatively large point spread function (PSF). 
The seven light curves from the discovery paper (jBakos et al.l 12003 ) were also obtained with the 
FLWO 1.2m and KeplerCam, and here we used a similar setup. 

Faulkes Telescope North (FTN): FTN is located on Mauna Haleakala in Maui, HI. We obtained 
four light curves from LCOGT's robotic 2.0 m telescope using the Spectral Instruments camera 
and a Pan-STARRS Z filter. The camera consists of a 4kx4k Fairchild Imaging CCD with a pixel 
scale of 0'.'304 pixel"^ (2x2 binning) and a FOV of 10f5 x 10^5. Exposure times ranged from 6 s 
to 10 s and a slight defocus was applied to the telescope in order to project the PSF onto a larger 
number of pixels, prevent saturation and increase the open shutter time relative to the overall cycle 
time. Three light curves were obtained during Season 3, and one during Season 4. 

Byrne Observatory at Sedgwick (BOS): BO^ is located at the Sedgwick Reserve near Santa 
Ynez, CA. We obtained two transit light curves of HAT-P-13b, both during Season 4, using the RC 
Optics 0.8 m remotely operated telescope at BOS. This telescope is equipped with a Santa Barbara 
Instrument Group (SBIG) STL-6303E camera containing a 3kx2k Kodak Enhanced KAF-6306E 
CCD with a pixel scale of 0'.'572 pixer^ (2x2 binning) and a 14!? x 9^8 FOV. We observed in the 
SDSS-i' filter, and exposure times ranged from 50 s to 80 s depending on atmospheric conditions. 
Due to the smaller aperture and short readout time (~10 s) at BOS, no defocusing was applied. 

At all observatories we gathered CCD images encompassing the target star IIAT-P-13. The 
moderately populated field surrounding the target provided several stars of similar brightness within 
the FOVs, to be used as comparison stars in the photometric processing. All data were reduced 
using standard routines for bias subtraction, dark current subtraction (when necessary), and flat- 
field correction. We extracted light curves with PyRAF using aperture photometry by dividing the 
flux of the target star by the weighted summed flux of several comparison stars in each image. Julian 
Dates of mid exposure we re recorded dur i ng th e observations, and later converted to BJD_TDB 
using the tools described in lEastman et~aD (|20inl i^. We optimized aperture sizes and the selection of 
comparison stars by minimizing the scatter of the resulting light curves, while iteratively removing 
5 a outliers. A total of 2 photometric outlier data points were removed from the collection of all 
light curves. All 10 new light curves are shown in Figures [1] and [21 and are listed in Table [TJ 

In addition to the 10 new light curves we obtained for this work, we also re-analyzed 12 light 
curves avai lable in the li t eratu re. Those include one light curve from Season 1 and six from Season 
2, all from iBakos et al.l (120091'). two l i ght cu rves presented by ISzabo et al.l (120101 ) from Season 3, 
and three light curves from lPal et al.l ( 2011 ) from Season 4. We adopted the previously published 
light curves as they were presented in their respective papers, and we only redid the fits in this 
work. Table [2] lists all 22 transit light curves included in our analysis. 



^Located at: 34.687604°, -120.039067°, 500m 



^Online tool for HJD_UTC to BJD_TDB conversion; |http://astroutils. astronomy.ohio-state.edu/time/hjd2bjd. html| 
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3. Analysis 
3.1. Analysis of all available data 

Our light curve fitting was done using the Transit Analysis Packag d fTAP: lGazak et alJboilh . 



Ford 


2005 


2006: 


Holman et al. 


2006; 


Collier Cameron et al. 


2007; 


Burke et al. 



20071 1. To account for possibl e temporally correlate d noise (e.g., iPont et al.ll2006l ) TAP uses the 



wavelet likelihood approach of lCarter &: WinnI (j2009l ). TAP has the ability to simultaneously fit 13 
parameters: orbital period (P), mid transit time (Tc), orbital inclination (i), orbital semi-major axis 
normalized by the host star's radius {a/Rg), planet to star radii ratio {Rp/Rg), two limb darkening 
coefficients {u\ and U2) for a quadratic limb darkening law, orbital eccentricity (e) and longitude of 
periastron (w). In addition, TAP fits a linear slope (5"), to account for a possible linear trend with 
time during the transit, a flux normalization factor (A), and two noise components: a temporally 
uncorrelated Gaussian "white" noise {a^j) and a time correlated "red" noise (0"^) (see equations 



32-34 of lCarter fc Winnll2009l ). where a power spectrum density of 1// is assumed. 

W e determined limb-darkening coefficient^ by interpolating over the grids of Claret ( 2000 . 



20041 ) and fixed these parameters in the analysis. Since e and uj are not well constrai ned by light 



curves alone Gaussian priors were assigned to these parameters using the values from IWinn et al 
mm -, e = 0.0133 ±0.0041, and oj = 210l^^ de grees. Therefore our model includes five parameters 



simultaneously fitted to all light curves (P, Tc, i, a/Rs, and Rp/Rg), and 22 sets of four parameters 
{S, N, aw, and ar) fitted individually to each light curve. We used jump rates of 25% for all free 
MCMC parameters, and ran 10 chains of 10^ steps each, discarding the first 10% of each chain 
before combining results of all chains. Each chain started from a different initial position 10 a away 
from the optimized parameter values. The best fit values, and upper and lower 1 a errors for each 
parameter were determined by taking the median, 15.9, and 85.1 percentile values respectively of 
the resulting a posteriori probability distri butions. In order to check the chains for non-conver gence, 
we calculated the Gelman- Rubin statistic ( Gelman et al. 2003 : Fordll2006 : Holman et al. 2006 ). The 
ratio of interchain variance to the intrachain variance was found to be within 10% of unity for each 
free parameter, giving no indication of non-convergence. Results of this analysis are shown in the 
bold row of Table [3l and we used those parameters for the over-plotted model in Figures [iHSl 

Table [2] includes parameters that indicate the quality of each light curve. The photometric 
noise rate (PNR) is defined as PNR=RMS/\/r, where the root mean square (RMS) is derived from 
the light curve residuals and F is the median number of cycles (including exposure time and any 
dead time such as readout time) per minute. Also listed are and ar, as fitted by the TAP. 



^ http://ifa.hawaii.edu/users/zgazak/IfA/TAP.html 



"111,^=0.5162 M2,v=0.2448, iii,fl=0.3971 it2,H=0.2977, mi, 7=0.2922 ?i2,/=0.3192, i^i,,. =0.3208 U2,i'=0.3124, 
Mi,z=0.2441 ?i2,z=0.3226 
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3.2. Seasonal Analyses 



We repeated this process treating the collection of all light curves from each of the three 
seasons 2-4 as separate subsets in order to look for possible variations in the system parameters 
from season to season. The resulting system parameters determined from all light curves and the 
three seasonal analyses can be found in Table [3l The bottom line of the table lists the parameters 
from lsakos et all ( 2009 ). for comparison. 



3.3. Refined Ephemeris 

We used the results from fitting all light curves in order to look for TTV and determine a 
refined ephemeris. For that end we analyzed each light curve separately by allowing only the mid 
transit time and the four light curve specific parameters (5, N, a^, and ar) to vary. The resulting 
mid-transit times for each transit event are listed in Table HI 

Once we determined the mid-transit times and the errors on those measurements, we then 
performed a linear least squares fit for a linear ephemeris, including P and a reference epoch Tcfi. 
Since we have some freedom in choosing the epoch for which Tc^o is fitted, we chose it to be during 
season 3, when the covariance between P and Tc^ is minimized, although we do not have a light 
curve of that specific transit event. The resulting parameters and their uncertainties are listed in 
the bold row of Table [3l We verified that the resulting cov(P, Tc^) is small enough and can be 
neglected when propagating the error bars to future (or past) mid transit times. 



Discussion 



The transit times O— C diagram showing the residuals from our linear ephemeris is presented 
in Figure IH and the residuals are listed in Table |4] as time difference, in seconds, and also after 
dividing by the mid transit times uncertainty, to show the significance of the difference. A close 
look at Figure m shows there i s only a singl e signi ficant outlier (12.7 min, 5.2 a), the earlier of the 
two transit events obtained bv lSzabo et al.l (j2010l ) during Season 3, on UT 2009 Nov 5. The linear 
fit to the mid-transit times produced = 45.57 with 20 degrees of freedom (DoF), and a reduced 
of Xred ~ 2-28. However, this value is highly affected by the ~5 a outlier from the UT 2009 
Nov 5 transit. If this single point is ignored we get = 19.36 with 19 DoF, and Xred ~ 1-02, but 
the difference in the resulting fit is small (Tc^o and P changed by —0.5 a and 0.1 a, respectively). 
The RMS of the 0-C residuals including the UT 2009 Nov 5 event is 211 s, and 144 s without 
including that event. 



Comparing our O— C diagram to the one presented by iPal et al.l ()201ll . see their Figure 2) 
shows a dramatic difference. Their figure shows that the mid transit times of the three light curves 
they obtained, during Season 4, strongly deviate from a linear ephemeris, by about 0.015 days, or 
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3- 18 a accordi n g to t he mid transit tii ne uncerta i nties they provide. The Hnear epher neris derived 



m 



Szabo et al.l (|20ld ) and adopted by iPal et al.l (j201l|) was b ased only on the s even iBakos et al 



( 200' 



^) transit times, from Seasons 1 and 2, and the two from Szabo et al. ( 2010 ). from Season 3 



Therefore the UT 2009 Nov 5 transit time from Season 3 heavily affected their derived epherneris. 

Here we have seven events from Season 3, including three observed by the FTN 2.0 m and 
two by the FLWO 1.2 m. Our larger number of observed transits from that season suggests that 
the UT 2009 Nov 5 mid transit time measurement is a single outlier, and that the data we have at 
hand are consistent with a linear ephemeris. 

We present the light curve from the UT 2009 Nov 5 event in Figure El over-plotted by our 
model from the analysis of all light curves and shifted to the best-fitting mid transit time for the 
UT 2009 Nov 5 event (left), and the mid transit time from the linear ephemeris (right). A close 
look at the light curve residuals, presented at the bottom part of both panels, shows it includes a 
few features, specifically during and before ingress, and during and after egress, and these features 
are clearly more pronounced in the right panel. It is possible that those features have affected the 
estimate of the mid transit time, and their origin could be astrophysical (although HAT-P-13 is 
not known to be an active star), or the result of correlated noise. 

We carefully examined the UT 2009 Nov 5 light curve, which consists of exposures alternating 
between V and R filters. We measured the mid transit time for the light curve observed in each 
filter independently, and found that they were both within 1.3 a of the mid transit time that we 
measured from the combined light curve, and close to 4 cr away from the linear ephemeris. We 
then applied a completely separate analysis of the combined light curve from both filters in which 
the parameters that determine the shape of the light curve (i, a/Rs, Rp/Rg, and Tc) were fitted 
and allowed to vary freely. The resulting fitted parameters obtained from the analysis of this single 
event are within 1.3 a from the values obtained from the Season 3 analysis. 



We ve rified that t he mi d transit times we derived here for the transits observed bv lBakos et al 



(|2009l ) and lPal et al.l (|201ll ) are consistent, within 0.7 cr, with the times derived by those authors. 
We measure a mid-trans i t tim e of rc,„=_i2=2455141. 552706 ± 0.001700, for the UT 2009 Nov 5 
transit from lSzabo et al.l (|20ld ) which is consistent with the value that they derive to within 0.3 a 



Our v alue of rc,„=25=2455249.447554 ± 0.001900 for the second transit obtained bv ISzabo et al 
|20ld ). from UT 2010 Feb 21, differs from the published value by 1.6 a. 



As already noted, the source of the large amplitude shift in the transit timing of the UT 2009 
Nov 5 event is unclear, and since our data do not include other transit events near that time it is 
difficult to rule out or confirm a physical process. If of astrophysical origin, it could be the result of 
an unusual physical process that affected also the shape of the light curve (see Figure [3]), and that 
is not seen during the other transit events analyzed here. HAT-P-13b transits observed by others 
close to that event, during October and November 2009 would be useful for shedding more light on 
this issue. 



We have attempted to look for a TTV signal in our O— C diagram, although it does not show 
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an excess scatter, besides the single outlier mentioned above. A parabolic fit to all 22 mid-transit 
times resulted in a value consistent with zero within 1 a for the quadratic coefficient, and with 
= 43.62 for 19 DoF, or x^ed = 2.30. Ignoring the UT 2009 Nov 5 outlier gives ^ = 19.31 for 
18 DoF, or x^ed = l-O^, and the RMS of the 0-C residuals are 210 s including the UT 2009 Nov 
5 event and 144 s without that event. Therefore, we could not identify a long term trend in the 
transit times. 

We also performed a p eriod analys i s on the resi duals from the linear ephemeris using the 
Lom b-Scargle (L-S) method (Lomblll976l : IScargldll982l ). looking for a possible low-mass perturber 



(e.g.. iHolman &: Murravll2005l ). The maximum peak of the periodogram was found to be similar to 



the maximum peak of periodograms in which the data were rearranged in a random order. More 
quantitatively, the strongest periodogram peak was at the 52nd percentile of a sample of strongest 
peaks in periodograms of 10^ random permutations, showing that no significant periodicity is seen 
in our mid t ransi t times O— C residuals. The large outlier from the UT 2009 Nov 5 event of 
Szabo et al.l (120101 ) was not included in this L-S analysis. 



S everal authors have presented predicted TTV behavior of HAT-P-13b transits (|Bakos et al 



20091 : iPavne &: Fordll201ll ). depending on the parameters of the second planet, HAT-P-13c. Our 
data put an upper limit on the ma ximum TTV aniplitu de of ~150 s during the 4 observational 
seasons. This reinforces the claim of lPayne &: FordI (|201ll ) that the eccentricity of the outer planet 
must be less than ~0.85, and the relative inclinations of the two planet's orbital planes must not 
be in the range 88°<ire/<92°. 

Of course it could be that there is yet another, short period low mass planet lurking in the 
system. To that end we note that a ~3 Earth mass planet orbiting at a coplanar orbit with twice the 
orbital perio d of HAT -P-13b will induce a TTV amplitude of 150 s (calculated using the methods 
presented in |Pai[201ol ). close to the detection threshold of our data. However, such a planet will 
induce also a 1 m s~^ radial velocity amplitude. Therefore, the sensitivity of our trans it timings to 



non transiting planets i n a 1:2 resonance is close to that of existing radial velocity data (iBakos et al 



200S 



Winn et al.ll2010l ). 



Our Season 2 light curves include data o nlv fromlBakos et al.l (|2009l ) . Comparing our results for 
that season with the parameters presented bv lBakos et al.l (|2009l ) shows they are in good agreement, 
and are consistent within 0.5 a (see Table[3]). This is an important validation of our analysis method 
using the TAP software. The uncertainties we derive are larger, though, by typically 40-60%, and 
330% for the orbital period. The latter can be explained by the additional value of the HATNet 
photometry in constraining th e period, bu t the former may indicate that our uncertainties are 
overestimated and/or those of iBakos et al.l ()2009l ) are underestimated. The error bars are also 
influenced by the fact that we did not include the UT 2008 Apr 25 event in the Season 2 analysis, 
since it was obtained during Season 1, although that light curve is partial and has a small impact 
on the parameters uncertainties. 



Thanks to our large amount of data, with a much longer time span, the errors on the system 
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light curve parame ters from the analysis of all 22 light curves are smaller than those obtained by 



Bakos et al 



((20091) by 15-35%, and more than a factor of 3 smaller for P. 



Our separate analysis of the data from each season (see Table [3]) shows a small shift, of ~ 2.5 a 
in i, a/Rs, and Rp/Rg, from Season 2 to Season 3, while the results for Seasons 3 and 4 are consistent 
within 1.2 a. The transit duration also shows a jump of ~ 2.8 a between Season 2 and Season 3, 
but the duration for Season 2, Season 4 and all 22 light curves is consistent to within 0.5 a. The 
shift in Season 3 may be influenced by the UT 2009 Nov 5 event, as the light curve (see Figure [3]) 
of this event may indicate a longer duration than other transits in our collection of data. The low 
significance of those shifts, and the fact that these parameters are correlated makes it difficult to 
draw any conclusion. If the UT 2009 Nov 5 light curve is excluded from the Season 3 analysis, 
then this jump becomes slightly less significant. The values of i, a/Rg, and Rp/Rs (i=81.67±0.87, 
a/i?s=5.29±0.31, i?p/i?s=0.0858±0.0023) if the UT 2009 Nov 5 event is excluded are consistent 
with the values of Season 3 when that event is included to within 0.4 a, and are less than 2 a away 
from Season 2. 



5. Summary 



We presented here an analysis of 22 HAT-P-13b transit light curves spanning four observational 
seasons, of which 10 were obt ained here and 1 2 were previously published. Contrary to the long 
term TTV signal suggested bv lPal et al.l (j201ll ) we find that the transit times are consistent with 
a linear ephemeris, while we identify a single transit time, from UT 2009 Nov 5, that significantly 
deviates from our linear model. The nature of this single deviation is unclear. The other light 
curve parameters do not show a large deviation compared to those fitted to light curves of the same 
season. Our large data set also allows us to refine the light curve parameters and transit ephemeris, 
which will be useful for future observational studies of this interesting system. 

This work demonstrates the use of a collaboration of ground-based 1 m class telescopes for 
transit timing monitoring, and that a large number of observations are required for thoroughly 
studying any TTV detection. This will und oubtedly be one of the goals of the future robotic 1 m 
class telescope networks, like LCOGT (e.g. IShporer et al.ll2010l ). Unlike expensive space missions 
as CoRoT and Kepler, small ground-based telescopes are easily accessible and their lifetime is 
not limited by the durations of space missions. Therefore, they will be an important resource in 
studying transiting planets orbiting bright stars in the decades to come. 



A.S. acknowledges support from NASA Grant Number NNX10AG02A. M.H. and J.W. grate- 
fully acknowledge support from the NASA Origins program through award NNX09AB33G. A. P. 
thanks the support of the ESA grant PECS 98073, and the Janos Bolyai Research Scholarship of 
the HAS. This paper uses observations obtained with facilities of the Las Cumbres Observatory 
Global Telescope. The Byrne Observatory at Sedgwick (BOS) is operated by the Las Cumbres 



Observatory Global Telescope Network and is located at the Sedgwick Reserve, a part of the Uni- 
versity of California Natural Reserve System. PyRAF is a product of the Space Telescope Science 
Institute, which is operated by AURA for NASA. 

Facilities: FTN (Spectral), FLW0:1.2m (KeplerCam) , LCOGT (BOS) 
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Table 1. Photometry of HAT-P-13b obtained in this work^ 



BJD.TDB Relative Flux 


Error 


Filter 


Telescope'' 


2455193.917958 


1.00015 


0.00089 


Z 


1 


2455196.836174 


0.99512 


0.00101 


Z 


1 


2455199.700075 


1.00361 


0.00166 


i' 


2 


2455231.842858 


1.00464 


0.00090 


z 


1 


2455275.603565 


1.00116 


0.00151 


i' 


2 


2455511.818399 


0.99610 


0.00163 


i' 


2 


2455613.891700 


1.00736 


0.00083 


z 


1 


2455616.787707 


1.00060 


0.00062 


i' 


3 


2455619.705125 


1.00093 


0.00061 


l' 


3 


2455622.690334 


0.99150 


0.00147 


i' 


2 



'^Only a sample is given here, the full table will be available in 
the online version of the manuscript 

^Telescope code is: 1 = FTN 2.0 m, 2 = FLWO 1.2 m, 3 = BOS 
0.8 m 



Table 2. Transit Observations of HAT-P-13b Analyzed in This Work 



Date'' 




Cycle Time 




a J- 


pnr'' 


Transit Part'^ 


Filter 


Telescope 


Reference 


(UT) 




(s) 


(%) 


(%) 


(% minute^"* 


(OIBEO) 


















Season 1 — 


2007/2008 








2008-04-26 


-204 


29 


0.19 


0.50 


0.14 


E O 


./ 


FLWO 1.2m 


Bakos et al. f2009'l 












Season 2 — 


2008/2009 








^iUUo-ii-UD 


-137 


29 


0.13 




0. 10 






r LjW \J 1 . zm 


JrSaKos ct al. ( ^iUUy ) 


2008-11-09 


-136 


29 


0.14 


0.56 


0.13 


OIBEO 


^' 


FLWO 1.2m 


Bakos ct al. (2009) 


2008-11-12 


-135 


29 


0.18 


0.13 


0.18 


B E O 


i' 


FLWO 1.2m 


Bakos ct al. (2009) 


2009-01-18 


-112 


29 


0.14 


0.14 


0.15 


OIBEO 


% 


FLWO 1.2m 


Bako^^L ^2009) 


2009-02-19 


-101 


29 


0.18 


0.55 


0.14 


O I B 


i' 


FLWO 1.2m 


Bakos et al. (2009) 


2009-05-09 


-74 


29 


0.15 


0.85 


0.14 


O I B 


i 


FLWO 1.2m 


Bakos et al. (2009) 












Season 3 — 


2009/2010 








2009-11-05 


-12 


132 


0.06 


0.57 


0.19 


OIBEO 


R&cV 


Konkoly 1.0m 


Szabo et al. (2010) 


2009-12-28 


6 


23 


0.17 


0.29 


0.11 


O I 


Z 


FTN 2.0m 


this work 


2009-12-31 


7 


22 


0.22 


0.58 


0.15 


O I B 


Z 


FTN 2.0m 


this work 


2010-01-03 


7 


44 


0.13 


0.19 


0.14 


O I B 


i 


FLWO 1.2m 


this work 


2010-02-04 


19 


30 


0.12 


0.46 


0.12 


OIBEO 


Z 


FTN 2.0m 


this work 


2010-02-21'* 


25 


411 


0.09 


0.58 


0.43 


OIBEO 


R 


Konkoly 0.6m 


Szabo et al. (2010) 


2010-03-20 


34 


34 


0.15 


0.38 


0.12 


O I B 


i 


FLWO 1.2m 


this work 












Season 4 — 


2010/2011 








2010-11-11 


115 


39 


0.16 


0.53 


0.14 


OIBEO 


i' 


FLWO 1.2m 


this work 


2010-12-27 


131 


39 


0.19 


0.75 


0.17 


OIBEO 


I 


Konkoly 0.6m 


Pal et al. (2011) 


2010-12-30 


132 


63 


0.08 


1.29 


0.17 


O I B 


R 


Konkoly 1.0m 


Pal et al. (2011) 


2011-01-28 


142 


28 


0.17 


1.44 


0.15 


O I B E 


R 


Konkoly 1.0m 


Pal et al. (2011) 


2011-02-21 


150 


31 


0.27 


0.55 


0.23 


O I B E 


Z 


FTN 2.0m 


this work 


2011-02-24 


151 


62 


0.10 


0.19 


0.11 


O I B 


i 


BOS 0.8m 


this work 


2011-02-27 


152 


94 


0.11 


0.37 


0.15 


OIBEO 


i' 


BOS 0.8m 


this work 


2011-03-02 


153 


34 


0.11 


0.65 


0.10 


B E O 


i' 


FLWO 1.2m 


this work 



^UT date at start of observation. 

^Photometric noise rate, calculated as RMS/v^, where RMS is the scatter in the light curve residuals and F is the median number of cycles 
(exposure time and dead time) per minute. 

OIBEO for Out-of-transit before ingress, Ingress, flat Bottom, Egress, and Out-of-transit after egress respectively. 

'^The long cycle time (and resulting PNR) is the result of ignoring the V band data of this event. Observations were o riginally taken whil e 
alternating between 'V and J? filters, but the V band data were plagued by large systematics and ignored in the analysis of lSzabo et all j2010l ) . 



REFERENCES 

Agol, E., Stcffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567 

Bakos, G. A., Howard, A. W., Noyes, R. W., Hartman, J., Torres, G., Kovacs, G., Fischer, D. A., 
Latham, D. W., Johnson, J. A., Marcy, G. W., Sasselov, D. D., Stefanik, R. P., Sipocz, B., 
Kovacs, G., Esquerdo, G. A., Pal, A., Lazar, J., Papp, I., & Sari, P. 2009, ApJ, 707, 446 

Burke, C. J., McCullough, P. R., Valenti, J. A., Jolms-Krull, C. M., Janes, K. A., Heasley, J. N., 
Summers, F. J., Stys, J. E., Bissingcr, R., FIccnor, M. L., Foote, C. N., Garci'a-Melendo, E., 
Gary, B. L., Howell, P. J., MaUia, F., Masi, G., Taylor, B., & Vanmunster, T. 2007, ApJ, 
671, 2115 

Carter, J. A., & Winn, J. N. 2009, ApJ, 704, 51 
Claret, A. 2000, A&A, 363, 1081 
— . 2004, A&A, 428, 1001 

Collier Cameron, A., Wilson, D. M., West, R. G., Hebb, L., Wang, X., Aigrain, S., Bouchy, F., 
Christian, D. J., Clarkson, W. I., Enoch, B., Esposito, M., Guenther, E., Haswell, C. A., 
Hebrard, G., Hellier, C, Horne, K., Irwin, J., Kane, S. R., Loeillet, B., Lister, T. A., Maxted, 
P., Mayor, M., Moutou, C, Parley, N., Pollacco, D., Pont, F., Qucloz, D., Ryans, R., Skillen, 
L, Street, R. A., Udry, S., & Wheatley, P. J. 2007, MNRAS, 380, 1230 

Eastman, J., Siverd, R., & Gaudi, B. S. 2010, PASP, 122, 935 

Ford, E. B. 2005, AJ, 129, 1706 

— . 2006, ApJ, 642, 505 

Gazak, J. Z., Johnson, J. A., Tonry, J., Eastman, J., Mann, A. W., & Agol, E. 2011, arXiv: 1102. 1036 

Gelman, A., Carlin, J. B., Stern, H. S., & Rubin, D. B. 2003, Bayesian Data Analysis, 2nd edn. 
(Chapman and Hall) 

Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288 

Holman, M. J., Winn, J. N., Latham, D. W., O'Donovan, F. T., Charbonneau, D., Bakos, G. A., 
Esquerdo, G. A., Hergenrother, C, Everett, M. E., & Pal, A. 2006, ApJ, 652, 1715 

Lomb, N. R. 1976, Ap&SS, 39, 447 

Mardling, R. A. 2010, MNRAS, 407, 1048 

Pal, A. 2010, MNRAS, 409, 975 



Table 3. Light Curve Parameters 





P 


i 


a/Rs 


Rpl Rs 


(BJD.TDB) 


(days) 


(deg) 







All 

Season 2 - 2008/2009 
Season 3 ~ 2009/2010 
Season 4 - 2010/2011 
Bakos et al. (2009) 



2455176.53880 ± 0.00034 

2454779.92895 ± 0.00072 

2455231.9464 ± 0.0012 
2455619.80708 ± 0.00085 
2454779.92979 ± 0.00038 



2.9162430 ± 0.0000030 

2.916305 ± 0.000033 
2.915952 ± 0.000063 
2.916203 ± 0.000048 
2.916260 ± 0.000010 



82.45 ± 0.46 

83-4 to.88 
81.16 ± 0.70 
82.05 ± 0.93 

83.4 ± 0.6 



5.52 ± 0.17 

5 86 

-0.34 

5.00 ± 0.23 
5.37 ± 0.36 
5.84 ± 0.26 



0.0855 ± 0.0011 

0.0838 ± 0.0019 
0.0882 ± 0.0020 
0.0857 ± 0.0021 
0.0844 ± 0.0013 
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Table 4. Mid-transit times of HAT-P-13b for the 22 light curves analyzed here 





T 




0-C 


0-C/o-T, 




(BJD.TDB) 


(s) 


(s) 




-204 


2454581.625065 


165 


-13 


-0.08 


-137 


2454777.012591 


106 


-79 


-0.74 


-136 


2454779.930096 


83 


30 


0.36 


-135 


2454782 843600 


197 


-206 


-1.05 


-112 


2454849.920092 


110 


44 


0.40 


-101 


2454882.000478 


204 


192 


0.94 


-74 


2454960.739683 


283 


248 


0.88 


-12 


2455141.552706 


146 


762 


5.22 


6 


2455194.035662 


198 


-51 


0.26 


7 


2455196.954496 


110 


173 


1.57 


8 


2455199.868674 


113 


-6 


-0.05 


19 


2455231.945420 


79 


-172 


-2.18 


25 


2455249.447554 


171 


232 


1.36 


34 


2455275.693121 


230 


178 


0.78 


115 


2455511.908538 


122 


155 


1.27 


131 


2455558.563744 


146 


-249 


-1.71 


132 


2455561.483633 


269 


66 


0.25 


142 


2455590.644415 


222 


-76 


-0.34 


150 


2455613.973903 


194 


-116 


-0.60 


151 


2455616.892899 


131 


122 


0.93 


152 


2455619.807862 


116 


11 


0.10 


153 


2455622.723514 


143 


-39 


-0.28 
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Season 3 - 2009/2010 




Fig. 1. — The five new light curves obtained in this work during Season 3, plotted in chronolog- 
ical order starting from the top, and offset in relative flux for clarity. The residuals after model 
subtraction appear below in the same order. The best fitting model, from the analysis of all light 
curves, is over-plotted in black. All dates are UT at the start of observation. 
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Season 4 - 2010/2011 




-0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 



Days from Mid Transit 
;. 2. — Same as Figure [H but for Season 4. 
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-0.05 0.00 0.05 

Days from Mid Transit 



-0.05 0.00 0.05 

Days from Mid Transit 



Fig. 3. — The UT 2009 Nov 5 light curve, which was obtained bv ISzabo et al.l (|20ld ). Our best- 
fitting model, obtained from TAP analysis of all 22 light curves is over-plotted in black. Left: The 
model is phased and shifted to the best-fitting mid-transit time for this event. Right: The model 
is phased and shifted to the best-fitting linear ephemeris from Table [3l 
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Fig. 4. — O— C plot including all 22 mid-transit time measurements. The black dashed lines indicate 
the 1 a errors on the predicted mid transit times by propagating the errors on P and from the 
linear fit. 



